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Abstract 

Inozemtsev models are classically integrable multi-particle dynamical systems re- 
lated to Calogero-Moser models. Because of the additional (rational models) or 
sin^ 2q (trigonometric models) potentials, their quantum versions are not exactly solv- 
able in contrast with Calogero-Moser models. We show that quantum Inozemtsev 
models can be deformed to be a widest class of partly solvable (or quasi-exactly solv- 
able) multi-particle dynamical systems. They posses AA-fold supersymmetry which is 
equivalent to quasi-exact solvability. A new method for identifying and solving quasi- 
exactly solvable systems, the method of pre-superpotential, is presented. 



1 Introduction 



In this paper we address the problem of the relationship among classical integrability, quan- 
tum integrability and quantum partial integrability (or quasi-exact solvability) within the 
multi-particle dynamical systems by taking a wide class of explicit examples, Inozemtsev 
models. We demonstrate that the Inozemtsev models (with degenerate potentials, that is 
non-elliptic potentials) can be made quantum partly solvable (or QES, quasi-exactly solv- 
able). For this purpose we present a simple new formulation of QES systems of single as well 
as multiple degrees of freedom. We also show that the notions of higher derivative (or non- 
linear or A/"- fold) supersymmetry and quasi-exact solvability are equivalent. In other words, 
Inozemtsev models provide plentiful examples of multi-particle A/'-fold supersymmetry. 
Inozemtsev models are generalisation of Calogero-Moser models 0, |TD 



rT|,|l^, |13l associated with the root system of BC type and A type. They belong to the 



category of 'twisted' Calogero-Moser models ||T2| , Like all the Calogero-Moser models 
they are classically integrable for all the four types of potentials, elliptic, hyperbolic, trigono- 
metric and rational in the sense that their equations of motion can be expressed in Lax pair 
forms. 

We start by asking a general and naturally vague question: to which extent does classical 
integrability imply quantum integrability in multi-particle dynamical systems? As is well- 
known, the converse, that quantum integrability always implies classical integrability, is 
true in multi-particle dynamical systems, since the quantum system is the h deformation 
of the classical one. We do not know a way to answer this question abstractly by starting 
from the pure notion of classical integrability, although some attempts have been made to 
construct quantum conserved quantities as a deformation of classical ones in the framework 
of perturbed conformal field theory (see, for example, |]15|). 

However, the experience of Calogero-Moser models, the widest class of integrable multi- 
particle systems ever known, tells that the classical integrability is very close to quantum 
integrability. The quantum integrability is proved universally, that is for all the root sys- 
tems including the non-crystallographic ones, for Calogero-Moser models with degenerate 
potentials |16|, [l^, namely those with trigonometric, hyperbolic and rational potentials. 

Thus we are naturally led to the question of quantum integrability of Inozemtsev models 
with degenerate potentials. Rational Inozemtsev models have an additional potential of sixth 
degree polynomial in q (coordinates), see ( |3.10|) , on top of the Calogero-Moser potentials. 



2 



Trigonometric Inozemtsev models have an additional sin^2g potential, see (|4.11|) , on top 
of the Calogero-Moser potentials. These additional potentials destroy the mechanism for 
providing quantum conserved quantities, the so-called 'sum-to-zero' condition of the second 
member of the Lax pair [|l^, 

The very interactions (g^, sin^2g, etc) that constitute obstructions for quantum integra- 
bility of Inozemtsev models are known to play essential roles in the quasi-exactly solvable 
|1T9| one particle quantum mechanics. This leads to a conjecture that at least a certain class 
of Inozemtsev models can be made quasi-exactly solvable. We demonstrate that supersym- 
metrisable (see ( p.5| ) and ( p.6| )) Inozemtsev models can be deformed to quasi-exactly solvable 
systems which are characterised by an integer deformation parameter Ai. It is also shown 
that the concepts of quasi-exact solvability and higher derivative or non-linear |^T| or 
M-fold [p2[ supersymmetry (with a typical relationship M = M. + 1) are equivalent. 



This paper is organised as follows. In section two we present the basic tool for investigat- 
ing quasi-exactly solvable systems which we call the method of pre-superpotential. We briefly 
summarise the classical Inozemtsev models in comparison with Calogero-Moser models. In 
section three we demonstrate the quasi-exact solvability of a single particle rational BC 
type Inozemtsev model based on a new method of employing a pre-superpotential W. This 
provides the most general single particle QES system with potential known to date. The 
equivalence of quasi-exact solvability and AT-fold supersymmetry is generally established. 



Other related notions, the "Bethe Ansatz" type equations p3[, ODE spectral equivalence 



2^ , |25| and Bender-Dunne polynomials are simply explained from the new point of 



view. Section four deals with the quasi-exact solvability of a single particle trigonometric 
BC type Inozemtsev model. In section five we discuss the rational A type Inozemtsev model 
with potential, which provides an example of a spontaneously broken A/'-fold supersym- 
metry. Through section six to section eight, various Inozemtsev models (rational BC type, 
trigonometric BC type and trigonometric A type) are shown to be quasi-exactly solvable and 
the generators of A/'-fold supersymmetries are identified. The QES of quantum Inozemtsev 
models is the consequence of exact solvability of quantum Calogero-Moser models and QES 
of the added single particle type interactions. The final section is devoted for comments and 
discussion. Appendix A gives the classical Lax pairs for the BC type and A type Inozemtsev 
models in the same notation as used in the main text. Appendix B is for the details of 
the lower triangularity of trigonometric Calogero-Moser interactions which are necessary for 
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establishing quasi-exact solvability of trigonometric Inozemtsev models. 



2 Basic tool and Classical Inozemtsev models 



2.1 Basic tool 

One basic tool for showing the existence of some exact eigenf unctions (quasi-exact solvability) 
is the following simple fact. Let W = W{q) be a real smooth function of the coordinate (s), 
then trivial differentiation formulas {p = —id/dq) 



2 W 

p e 



+ 



-W 



E2 W 
Pje 



■E 



fdW' 



\dqj 



+ 



dqj 



imply that is an eigenfunction of the Hamiltonian H with eigenvalue 0: 



He 



w 



0, 



H = -p' + - 
2^ 2 



+ 



dq^ 



-I r 1 ^ 



+ 



dq] 



(2.1) 



so long as it is square integrable 



1 e'^rf^g < oo. (2.2) 

This is the simplest example of quasi-exact solvability. Looked differently, one might say 
this is a property of 'factorised' Hamiltonians 

.dW\ I 



H 



dW\ ( dW 
P + i- 



H 



dq 2 f^, V % 

together with a differential operator(s) that annihilate the state 



Pi 



. dW 

' 9qi , 



(2.3) 



p + i- 



0, 



Pi + i- 



0, J = l, 



(2.4) 



\ '^q ) ' V ^qj J 

This fact can also be considered as the very base of supersymmetric quantum mechanics. 

This gives the ground state wavefunction of the quantum Calogero-Moser models. In 
other words, all the (quantum integrable) Calogero-Moser models can be described by pre- 
superpotentials W 0, |T8|. To sum up, if a Hamiltonian can be expressed in terms of W 
as ( |2.1|) or (|2.3| ) up to a constant, the existence of one eigenfunction is guaranteed save the 
square normalisability. Throughout this paper we call function W a pre-superpotential . 
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2.2 Classical models 



Here we present classical Inozemtsev models together with Calogero-Moser models for com- 
parison. The dynamical variables are canonical coordinates {qj} and their canonical conju- 
gate momenta {pj}- We denote them by r-dimensional vectors q and p with standard inner 
product: 

g = (gi, ■ ■ ■ G R^ p = {pi, ■ ■ ■ ,Pr) e W , p^ = p-p = pl^ V p^.. 

As is well-known for a root system A (rank r) and four types of potentials, elliptic, trigono- 
metric, hyperbolic and rational, Calogero-Moser and Inozemtsev models are classically com- 
pletely integrable. In this paper we discuss only those models based on classical root systems, 
that is A type and BC (and D) type and degenerate potentials, that is trigonometric, hy- 
perbolic and rational potentials. Since algebraic structures are almost the same for the 
trigonometric and hyperbolic potential models, we discuss trigonometric case as a repre- 
sentative. Among various types of Inozemtsev models we focus our attention to the 
supersymmetrisahle models, namely to those models whose Hamiltonians can be collectively 
expressed in terms of a pre-superpotential W = W{q) ||TH| as 

:V 

or 'factorisable' at the classical level: 



1 ^ / ,dW\ ( .dW\ 

Each specific model in this class is given by a choice of W , which are listed as below. 

2.3 BC type Calogero-Moser models 

The rational model pre-superpotential W 

r r 

Wc-M = Qm {log l^i - + log l^i + (lk\] + 9sY. log l^il' (2-7) 

j<k j=l 

contains two real coupling constants Qm for the middle roots (length^ =2) and gs for the 
short roots (length^=l). The trigonometric model pre-superpotential W 

r 

Wc-M = 9m J2 {log I sin(gi - g/c) I + log I sin(gj + qk)\} 
j<k 

r r 

+9s log I sin qj \+ gij^^^sl sin 2qj \ , (2.8) 



has one more coupling constant than the rational case, Ql for the long roots (length^=4). 
For the rational potential, the long roots and short roots are essentially the same. When 
Ql = and (75 = the models belong to the Dr root system, li = and gs ^ { 9l ^ 
and gs = 0) the models belong to the Br (Cr) root system. Throughout this paper we put 
the scale factor in the trigonometric functions to unity for simplicity. 

2.4 A type Calogero-Moser models 

For A type models, it is customary to express the roots by embedding in a space with one 
higher dimensions. We will discuss A^-i models with r degree of freedom. Since all the roots 
have the same length, A type models have only one real coupling constant g. The rational 
model pre-superpotential W is given by 

r 

Wc-M = gJ2^^^ h - (2.9) 

j<k 

whereas the trigonometric pre-superpotential W reads 

r 

M/'c-M = ^^log|sin(gj -gfc)|. (2.10) 

j<k 

2.5 BC type Inozemtsev models 

The rational supersymmetrisable BC type Inozemtsev model 0, § has two more real coupling 
constants, a and b, than the corresponding Calogero-Moser model, 

W = -j: + ^qfj + Wc-m(2.7), (2.11) 

which leads to degree six polynomial potentials. The trigonometric supersymmetrisable BC 
type Inozemtsev model [|, ^ has also two more real coupling constants, a and b, than the 
corresponding Calogero-Moser model, 

^ = E (-^cos2g, + ^log|cotg,|^ +Wc-m{2.8). (2.12) 

2.6 A type Inozemtsev models 

The rational supersymmetrisable A type Inozemtsev model has two more real coupling 

constants, a and b, than the corresponding Calogero-Moser model, 

^ = E (1^1 + ^^') + W^c-m(2.9), (2.13) 
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which leads to degree four polynomial potentials. The trigonometric A type Inozemtsev 
model has only one more real coupling constant, a , than the corresponding Calogero- 

Moser model, 

^ = E (-|cos2g,) + iyc-M(2.10). (2.14) 

It should be emphasised that these additional interactions of the supersymmetrisable 
Inozemtsev models are all 'single particle' type. In the following three sections §|^§| we 
will investigate the characteristic single particle dynamics of supersymmetrisable Inozemtsev 
models. 

3 Rational BC type Inozemtsev model with one degree 
of freedom 

Let us start with a pre-superpotential W = W{q) which is decomposed into two parts: 

W = Wo + Wu (3.1) 
'a 4 6 
4^ +2 



^0 = -( V + ^g' 1+^5 log |g|, a>0, gs>0, (3.2) 



M 

W, = J2log\q'-U (3.3) 

k=l 

in which A4 is an arbitrary non-negative integer and {^/tj's are distinct but as yet undeter- 
mined parameters. The first part Wq corresponds to the 'single particle interactions' of the 
BC type rational Inozemtsev model, which is an even function of q. The added part gives 
rise to an arbitrary polynomial in of degree Ai in = Ylk^iiQ"^ — ^k)- Since 

\ dq J dq^ \ dq ) dq^ dq dq \ dq ) d(f 

we will evaluate the terms containing W\. 

^diW^dm^ ^ (diWx^ ^ d'^Wi 



dq dq \ dq I dq^ 



r , „ , ^5! ^ 2g ^2 2v- 1 1 

2 -{aq^ + hq) + ^-^\Y: + E + 8^ E ^ 

I q It^^q"- ik fctl <f-ik t^i<f- ik - 



(3.4) 



This is a meromorphic function in q^ with at most simple poles. We demand that the residues 



of the simple poles, q = ^k-, k = 1, . . . , M. should all vanish |2^, which results in a set of 
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rational ("Bethe ansatz" type) equations for {^fc}'s: 

2 {-{ae^ + b^k) +9s} + l + ^^k E = 0> k = l,...,M. (3.5) 

i^k '5' 

Then expression (|3.4| ) becomes a linear polynomial in g^, which is easy to evaluate 

M 

(0) = -4aMq^ - 4bM - 4a ^ ^k- 

k=l 

Thus we arrive at 

i^j +^ = i^j +^-^-M,^-^E,, (3.6) 

El = 2bM + 2aJ2^k- (3.7) 

fc=i 

It should be emphasised that except for the constant term Ei the expression ( |3.6| ) is inde- 
pendent of the parameters {^/cj's introduced in Wi. 

This means that, for each set of solutions {■Cfcj's (with real J2Ck and up to the ordering), 
we have an eigenstate 

M M 

= e^o n i<l" - = q'' e'^t'^^+l'^^) J] {q' - Ck) (3.8) 

k=l k=l 



with eigenvalue Ei, (p.7|), of the Hamiltonian 



2^ 2 



dq I dq^ 



2aMq^, (3.9) 



-p + -q [aq + b) + — a(- + 2M + gs)q - - 



It has g^, g^, and 1/g^ potentials and a part of the coefficients of the quadratic potential 
is quantised. Because of the singular centrifugal term 1/g^, we restrict the function space to 
the half line, (0, +oo). The restriction on the coupling constants, a > is for securing the 
square integrability of at q = +oo and gs > for finiteness at g = 0. 

The above result implies that by adding a single term — 2aA^g^ to the Hamiltonian the 
single particle BC type rational Inozemtsev model can be made quasi- exactly solvable, that 
is a finite number of eigenstates together with their eigenvalues can be obtained exactly 
by algebraic means. The very term — fg^ in Wq that obstructs quantum integrability is 
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instrumental for the introduction of the additional term —2aM.q^. The eigenfunction 
of the above quasi-exactly solvable system belongs to a "polynomial space" 



^2k 2M 
,...,(/ 



Vm = Span ^l,q ,.. .,q 
In other words, the Hamiltonian ( |3.10D leaves this polynomial space invariant: 



(3.11) 



(3.12) 



Therefore these "polynomial space" can be called the exactly solvable sector of the system. 
It is elementary to see that the "polynomial space" Vm is annihilated by an A/" = + 1 st 
order differential operator P^: 



PnVm 



N-l ( 

n 

k=0 \ 

0, 



D 



t{Ar-iy 



D 



D = p + i 



dWp 
dq 



(3.13) 
(3.14) 



Since Pj^ is an A/" = + 1 st order differential operator, it is obvious that Vm gives the 
entire solution space of a differential equation 

Puy = 0. 

This differential operator together with its hermitian conjugate defines a higher derivative 
| 20| or non-linear or an Af-fold p2[ supersymmetry generated by 



Q = P^^\ Q^ = pU, 



(3.15) 



in which ip and ip'^ are fermion annihilation and creation operators. The "polynomial space" 
Vm is characterised as the zero-modes of Q and 



QVm = Q^Vm = 0, 



(3.16) 



which is the generalisation of the property of the ground state of the ordinary (A/" = 1) 
supersymmetric quantum mechanics. The second equality Q'^Vm = is trivial since Vm has 
zero fermion number. 

The structure of the exactly solvable sector can be better understood by making a simi- 
larity transformation of H by e^'^ (see for example): 

H 



2 dq dq 

+ {a(f + hq-—)^- 2aMq^. 
2 q dq 



(3.17) 



Then the above "polynomial space" is mapped to a genuine polynomial space 



Vm = Span 
whose invariance under H ( p.l7 ) 



,2M 



(3.18) 



is rather elementary to verify. If one substitutes an expansion 



M 
k=0 



2k 



"0 



into the eigenvalue equation 



H<il = Em 

one obtains a three term recursion relation for {a^j's 

{k + 1)(2A; + 1 + gs)ak+i = {2kh - E)ak + 2a{k -M- l)ak-i. (3.19) 
This determines as a polynomial in E of degree k which is a Bender-Dunne polynomial 



26| in the naivest sense. The condition a^+i = gives the characteristic equation of H: 

aM+i = ^ {2Mb - E)aM - 2aaM-i = ^ det{H - E) = 0. (3.20) 



In the rest of this section let us discuss the relationship between quasi-exact solvability 
and A/"- fold supersymmetry in the general context. This applies to the other cases discussed 
in the following sections as well. The exactly solvable sector of a quasi-exactly solvable theory 
is characterised by its "polynomial space" 



Vm = Span 

which is invariant under Hamiltonian 



l,/i, 



gen 



HgenVM ^ Vm- 



(3.21) 



(3.22) 



In these formulas the subscript "gen" in H and W stands for 'generic' and the function 
h = h{q) need not be a polynomial in q. (For example, in the trigonometric BC type 
Inozemtsev model (see section ^) h{q) = sin^g.) It is straightforward to verify that the 
"polynomial space" Vm is annihilated by an A/" = + 1 st order differential operator P^: 



Pm = Y[iDgen + ikE{q)), Dgen=P + i 

k=0 

PmVm = 0. 



gen 



(3.24) 
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As above, the "polynomial space" Vm gives the entire solution space of the differential 
equation P^y = 0. One could summarise the situation as the exactly solvable sector of 
a quasi- exactly solvable dynamics is characterised as the states annihilated ( |3.16|) by the 
generators Q and ( |3.15|) of an M-fold supersymmetry. 

On the other hand, let us suppose that one has a pair of Hamiltonians Hg^n and H^^^ 
which are intertwined by Pj\f p2| , p4| , p5| : 



Pj^Hgen - H+^Pm = 0. (3.25) 

Let Vm be the space of solutions of the differential equation P_\fy = 0, which is finite 
dimensional. Then from (|3.25|) one obtains Pj^H^s^Vm = and thus deduces that the finite 
dimensional space Vm is invariant under Hgen- H^chVm ^ '^M, ( P-22|) . One could summarise 
this as the quasi-exact solvability of -ffgen is a consequence of the M-fold supersymmetry and 
the intertwining relation ( p.25| ). The spectral equivalence of -ffgen and i^g^n holds outside of 



Vm as in the ordinary (TV = 1) supersymmetric quantum mechanics. 

4 Trigonometric BC type Inozemtsev model with one 
degree of freedom 

Here we consider a one dimensional quantum mechanical system with the following super 
potential in a finite interval [0,7r/2]: 

W = Wo + Wi, (4.1) 
d b 

Wo = --cos2g + - log I cot g| + gf^log | sing|, (4.2) 

9s >0, gs>^>0, (4.3) 

M 

W, = Y.log\sm'q-U (4.4) 

k=l 

in which Wq part is obtained by retaining the single particle part of the trigonometric BC 



type Inozemtsev model ( |2.12D . It is an even function of q and it reduces to a well-known 



"double sine-Gordon" quantum mechanics if only the first term |cos2g is kept. Here 



we have not included the long root term g^ log | sin 2q\ in (|2.8| ) since it can be expressed as a 
linear combination of log | cot q\ and log | sing| terms. As in the previous section, we evaluate 
the terms containing Wi in {dW/dq^ + d'^W/dq^: 

b \ ^ 2 sin g cos q 



2 a sin 2g r 



\ M 

+ 9s cot q\Y^ 

/ k=l 



sin2g ' '"'"""^It^. sin^ g - 
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+8 sin^ q cos^ — 



1 



M 



+ 2 COS 2g ^ — 



which is a meromorphic function in x = sin^ q with at most simple poles: 
(4.5) = 2(Aax{l-x)-h + 2gs{l-x))Y.^- 

k=l ^ 



(4.5) 



M 



1 



+ 2(l-2x)^ _ . 

k=i ^ 



(4.6) 



As in the previous case we demand that the residues at the simple poles x = ^fc, = 1, . . . , 
should all vanish. This requires that {^fc}'s should obey a set of rational equations 



i^k ^' 



0, k = l,...,M. (4.7) 



Then expression (|4.6|) becomes a linear function in x = sin^ q which is easy to evaluate 



M 



) = -8aM sin^ g _ 8a ^ ^fc - 4:M{gs + M). 

k=l 



Thus we arrive at 



'dW' 



+ 



dq I dq 



'dWo 
, dq 

M 



dq^ 



-8aM sin2g-2Ei, 



E, = AaJ2^k + 2M{gs + M). 

k=l 



(4.8) 
(4.9) 



Again, except for the constant term the expression ( |4.8| ) is independent of the parameters 
{^fcj's. Thus for each real solution of ( ^.7| ), we have an eigenfunction 



M 



M 



e- = l[{sm^ q-^k) = {sin qys {cot q) 2 e-^''"'^'^l[{sm^q-Ck 

k=l k=l 



(4.10) 



with eigenvalue Ei, ( [4.9|) , of the Hamiltonian 

'dWo\ 



H 



1 2 1 

2^ 2 



dq 



+ 



dq^ 



AaM sin^ q. 



(4.11) 



The eigenfunction ( [l.iup is square integrable in [0, 7r/2] due to the restriction on the pa- 
rameters (|4.3| ). In other words, the above Hamiltonian ( [4.11| ) is quasi-exactly solvable. The 
eigenfunction ([4.101) belongs to a "polynomial space" 



Vm = Span 1, sin g, . . . , (sin q) , . . . , (sin q) 



\2M 



(4.12) 
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which is invariant under the action of the Hamiltonian HVm ^ "^m- K is easy to see that 
Vm is annihilated by an A/" = + 1 st order differential operator Pj^: 

Px = Y[{D + i2kcot2q), D=p + t—^, (4.13) 

fc=0 

P^Vm = 0. (4.14) 

Thus the general statements in the previous section concerning the quasi-exact solvability 
and A/'-fold supersymmetry also hold in this case. 

In order to investigate the structure of the exactly solvable sector we make as before the 
similarity transformation of H by e^°: 

H = e-^»i/e^° = V-?^|--4aA<sin2g, 

2 oq oq 

= -p — (asin2g \-gs cot q)— 4aAlsin q. (4.15) 

2 sm2g oq 

The eigenfunction of H in the polynomial space Vm = Span 1, sin^ g, . . . , (sin g)^-'^ , can 
be obtained by substituting an expansion \E' = Y.k^oOik{sinqy'' , ao = 1 into the eigenvalue 
equation if\E' = E^^. This again leads to a three term recursion relation for "Bender-Dunne" 
polynomials {afc(-E)}'s: 

{k + l){2k + l-b + 2gs)ak+i = (2A;(1 - 2a + 2gs) - E) ak + 4a(A: -M- l)ak-i. (4.16) 

Again ak{E) is a polynomial in E of degree k and the condition om+i = gives the charac- 
teristic equation for H: 

c^M+i = 0^ {2M{1 - 2a + 2gs) - E) - ^auM-i = det(i^' - E) = Q. 

5 Rational A type Inozemtsev model with one degree 
of freedom 

This is an interesting example which fails to achieve quasi-exact solvability due to the lack of 
square integrability of the eigenfunction. It is interesting to know how far the algebraic pro- 
cedures go in parallel with the previous cases. We start with the following pre-superpotential 
W which is obtained by retaining the single particle part of W in ( |2.13| ): 

n h ^ 
W = Wo + Wu Wo = -q' + -q\ W^i = J] log |g - ^1- (5.1) 

^ ^ k=\ 

13 



This is cubic in q and leads to a quartic potential of g, see ( p.5|) . The terms containing Wi 
in {dW/dqY + d'^W/dq^ are: 



k=i 9 ~ fc<« i-ikQ-ii 



(5.2) 



From the requirement of vanishing residue at g = ^fc, we obtain rational equations 

1 



i^k 



0, k,l = l,...,M, 



(5.3) 



and the expression ( |5.2|) reads 



3) = 2aMq + 2aY,^k + 2bM. 



k=l 



Thus for each real solution {C,k} of ( p.3| ) we obtain an "eigenfunction" 



A;=l 



M 

n(9-e.) 

fc=l 



(5.4) 



of a Hamiltonian 



i7 



-p + - 
2^ 2 



. dq 



+ 



dq^ 



'-q^{aq + 6)^ + a(A< 



+ aMq, 



2^^+ 2' 



with energy 



M 



E 



k=l 



(5.5) 



(5.6) 



Let us recall the simple facts about the limiting case of a = and b = —u, u > 0. Then the 
Hamiltonian ( ^.5| ) becomes that of the simple harmonic oscillator with angular frequency u 
and the equations ( ^.3|) determine {^fe}'s as the zeros of Hermite polynomials p8| , with scaling 
by ^/uJ. This results in the well-known eigenfunction with Hermite polynomials ( ^.4|) and 
the spectrum E = ujA4 ( |5.6| ). (The zero-point energy (X'/2 is contained in the Hamiltonian 
(|5.5|).) Thus at least for 6 < and \a/b\ ^ 1, that is the quartic and the accompanying cubic 
terms in the potential can be considered as 'perturbations', it is expected that the equations 
( |5.3| ) have real solutions and the above solution generating method would work. However, 
the "eigenfunction" ( ^.4|) is not square integrable in the region (— oo, +oo) for whichever 
choice of the sign of a 7^ 0. One might be tempted to restrict the region to a half line, say 
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(0, +oo), by introducing a singular potential at the origin, for example, by adding a term 
(7log|g| to Wq, as in the example in section However, this cannot remedy the situation 
because of the wrong parity of the aq^ term. 

The J\f = Ai + 1 st order differential operator annihilating the "polynomial space" 



V 



M 



Span 



M 



has a very simple form: 



Pm = D^, D=p + i- 



dq 



(5.7) 



In this case the AT-fold supersymmetry generated hy Q = Pj\fip^ and = Plfi^ is sponta- 
neously broken for a 7^ 0. 

The difference in the algebraic structure from the quasi-exactly solvable case discussed 
in section § becomes clearer by making the similarity transformation of H by : 



H = e-'^^He^'' = lp'~^^ + aMq, 

2 oq oq 

1 d 

= -p^ - iaq^ + hq)— + aMq. 

2 oq 



(5.9) 



This maps q^ to q^^^, q^ and q^ ^: 



Hq' = a{M - k)q'+' - bkq" - UMl-^qk-^ 

Since the last term is g'^"^ instead of q''~^, the three term recursion relations for the co- 
efficients {afcj's in a series solution = Y.k^oCekQ'' , chq = 1 for the eigenvalue equation 
= E^^ do not hold any more. The {ak{E)}^s are no longer degree k polynomial in E. 

We will not discuss "trigonometric A type Inozemtsev model with one degree of free- 
dom", since the single particle part of ( p.l4|) is simply W = — |cos2g, that is the "double 
sine-Gordon" quantum mechanics. This is a well-known example of quasi-exactly solvable 
dynamics and is a special case of the model treated in section ^. 



6 Rational BC type Inozemtsev model 

Following the results of the single particle case in section ^, we consider the following Hamil- 
tonian 



^ i=i ^ 



dq. 



+ 



dq] 



- 2aM y: q] 



(6.1) 
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in which Ai is an arbitrary non-negative integer and Wo is given by (|2.11|) : 

Wo = - E f + ^^l) + ^c-M, (6.2) 



Wc-M = 9Mj2{^^s\(lj-(lk\+'^og\qj + qk\} + gsY.^og\qj\, 

j<k j=l 

a>0, gs>0, gM>0. (6.3) 

The Hamiltonian as well as Wq are Coxeter (Weyl) invariant of the BCr root system. 
The only difference with the classical supersymmetrisable Inozemtsev model is the added 
quadratic terms proportional to Ai. We will consider the model in the principal Weyl cham- 
ber 

gi > > ■ ■ ■ > gr > 0. (6.4) 

A special case of this Hamiltonian with one free parameter other than Ai (plus an invisible 
overall scale factor) was discussed in [p^ . 

In order to show the quasi-exact solvability of the Hamiltonian (|6.1| ) we have to demon- 
strate that a certain "exactly solvable sector" is invariant under H. As before let us first 
define a "polynomial space" 

Vm = Spano<„^<^, ,<^.<, [{qlr " " " (q-P " " " iq^r] e^°. (6.5) 

The "exactly solvable sector" is the permutation {qj <-> g^) invariant subspace of Vm- 

V^^^ = {veVM\9v = v, \/geGBc}, (6.6) 

in which Gbc is the Coxeter (Weyl) group of the BC root system. This fact can be seen 
easily, as in the single particle case, by similarity transformation 

1 r r QTi/ Q r 

2fr{ ' dqj dqj fr{ ' 

2,tl ' p{ % dqj ^tl ' % p{ ' 

Vm = Spano<„^.<^, ,<^.<,[(g?)"^---(g|)"---(g2)"^]. (6.8) 

The proof of the invariance of Vm^^' under the added terms Yl^j=i{aqj +bqj)^ — 2aAi Z]j=i Qj 
is essentially the same as in the single particle case. As for the Calogero-Moser part, W^c-m, it 
always decreases the power I]j=i'^i by one unit. The Coxeter (Weyl) invariance is necessary 
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and sufficient so that the result remains a polynomial, without developing unwanted poles. 
Thus the invariance of V^'^ under H and the quasi-exact integrability of H is proved. It is 
obvious that the elements of are square integrable since 

j=i j<:k 

and the restriction on the parameters ( |6.3|) and the integration region ( |6.4|) secure finiteness 
at infinity and at the boundaries of the Weyl chambers. 

Thus we have shown that the rational BC type Inozemtsev model can be made quasi- 
exactly solvable by adding properly quantised quadratic terms —2aM. Z]5=i 1j- The above 
"polynomial space" Vm is annihilated by the following r different JV = Ai + 1 st order 



commuting differential operators Pjip'- 










= n(A + 


.k 




(6.10) 








= 0, j,k = l,...,r. 


(6.11) 


The A/'-fold supersymmetry is 


generated by 








Q 




Q'- 


- tipi^^n^ 


(6.12) 



j=i i=i 
in which ijjj and ijjj are the annihilation and creation operators of the j-th fermion. 

7 Trigonometric BC type Inozemtsev model 

Following the results of the single particle case in section ^, we consider the following Hamil- 
tonian 



^ i=i ^ j=i 



+ 



2 



4aA^^sin2gj., (7.1) 



in which Ai is an arbitrary non-negative integer and Wo is given by ( |2.12| ): 

^0 = E (-^cos2g, + ^log|cotg,| I +iyc-M, (7.2) 



Wc-M = 9m J2 {log I sin(gi - gfc) I + log I sin(gj + qk)\} 
j<k 

r 

+5-5 log I sin I- (7.3) 
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All the parameters a, b, qm and gs are real and they satisfy 



^5 > 0, C/Af > 9s > ^> 0. 



(7.4) 



Let us recall here the argument in section ^for dropping ql log | sin2gj| term in favour of 
h log I cot qj I term. The Hamiltonian as well as Wq are Coxeter (Weyl) invariant of the BCr 
root system. The only difference with the classical supersymmetrisable Inozemtsev model is 
the added sin^ qj terms proportional to M.. We will consider the quantum mechanical model 
in the principal Weyl alcove 

7r/2 > gi > g2 > ■ ■ ■ > > 0, (7.5) 

due to the periodicity and Coxeter (Weyl) invariance of the model. As in the single particle 
case, we will demonstrate that a certain "exactly solvable sector" is invariant under H. As 
before let us first define a "polynomial space" 



V^ 



[sin^ gi)"i . . . (sin^ g,)"^ ■ ■ ■ (sin^ q,y- 



(7.6) 



cos2(^n^.gj) 



(7.7) 



— ^P^^O<nj<M, l<j<r 

or equivalently: 

Vm = ^P^^-M<n'.<M, l<j<r 

The "exactly solvable sector" is the permutation {qj ^ q^) invariant subspace of Vm'- 

V^^'' = {veVM\gv = v, WgeGsc}, (7.8) 

in which Gbc is the Coxeter (Weyl) group of the BC root system. All these functions are 
square integrable, since the integration region is finite and the possible singularities in e^" 
at the boundary ( |7.5| ) 



(7.9) 



^ _ e tE;=i-«2,, J|(sing.)95(cotg,)t Y[[sm{qk - qi) sin(gfc + qi)]^^^ , 

3=1 k<l 

are taken care of by the restriction on the parameters ( [7.41 ). The quasi-exact solvability can 

be shown by the similarity transformation 

1^2 ^dWo d 
- > n — y — 

d 



H 



nl^Pj ^ 2^(asm2gj 

^ ' 1 dqj dq, ^tl 



h 



2 ^ 
^ i=i 3 



V 



'^V^'^0<nj<M, l<j<r 
^P^^-M<n'.<M, l<j<r 



[sin^ gi)"i ■ ■ ■ (sin^ qj)""^ ■ ■ ■ (sin^ qr 



cos2{J2n'jqj) 

i=i 



(7.11) 
(7.12) 
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The corresponding "exactly solvable sector" is the permutation {qj ^ qu) invariant subspace 
of Vm- 



m\9v = V, 



(7.13) 



As for the Calogero-Moser part, the above Hamiltonian is lower triangular [|T6| in the basis 
( B.6 ) of V'^^ . The lower triangularity is stronger than the invariance of the polynomial 



space under the Hamiltonian. An outline of the proof is given in Appendix B. As for the 
added part 

-t (asin2g, |- - 4aA^ E sin^ , 

the proof of the invariance of the polynomial space (|7. 1 1| ) is essentially the same as in 
the single particle case. Thus the quasi-exact solvability of the trigonometric BC type 
Inozemtsev model is established. 

The above "polynomial space" Vm ( |7.6| ) is annihilated by the following r different TV = 



+ 1 st order commuting differential operators P 



U). 

Af ■ 



P. 



dW 

\{{Dj + 2ik cot 2qj ) , Dj= pj + i , 
0, [P^\plf)]=0, j,k = l,...,r. 



(7.14) 
(7.15) 



The A^-fold supersymmetry is generated by 



i=i i=i 
in which ipj and ip] are the annihilation and creation operators of the j-th fermion. 



8 Trigonometric A type Inozemtsev model 



This has a much simpler Hamiltonian than the previous one: 



+ 



d^Wo 
dq] 



AaM ^ sin^ qj, 



in which Ai is an arbitrary non-negative integer and Wq is given by (|2.14|) : 



^0 = -^Ecos2gj + H^c-M, 
^ i=i 



Wc 



C-M 



^ E log I si 



sm 



Qk) 



j<k 



(8.2) 
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There are only two real coupling constants a and g and we require 

9 > (8.3) 
for square integrability of the eigenfunctions of the form (|8.5| ). The Hamiltonian and Wq are 



invariant under any transpositions {qj,Pj) ^ {qk,Pk), which form the Coxeter group of Ar-i. 
Thus we consider the quantum mechanics in the fundamental Weyl alcove: 



TT > gi > g2 > ■ ■ ■ > Q'r > 0. 



Reflecting the simple form of the Hamiltonian, the "exactly solvable sector" which is left 
invariant under H has a simpler structure than those in the previous cases. Let us first 
define a space of truncated Fourier series with two units: 



V 



M 



Span_ 



M<nj<M, l<j<r 



The "exactly solvable sector" is the permutation [qj ^ q^) invariant subspace of V 

V^^ = {veVM\gv = v, WqeGa}, 



.5) 



(8.6) 



in which Ga is the Coxeter (Weyl) group of the A type root system. The quasi-exact 
solvability of the Hamiltonian ( ^.1|) is again easily verified by the similarity transformation: 



H 



p-Worr Wo _ ^ V „2 9Wo d 



AaM sin^ g^, 



g dq, dq, 



V 



M 



Span 



— a sm iqj — 4a Al 2^ sm , 

oqj 



1.7) 



-M<nj<M, l<j<r 



The lower triangularity of the Calogero-Moser part of the Hamiltonian in the basis ( [B.l|) of 

(|8.8| ) was proven originally by Sutherland 0. For the additional part 

^ Q ^ 

—a ^ sin 2gj ^ AaAi ^ sin^g^, 

j=i j=i 

the proo f that it leaves the space of the truncated Fourier series ( |8.8|) invariant is rather 

elementary. 

The above "polynomial space" Vm ( |8.5| ) is annihilated by the following r different Af = 
2Ai + 1 st order commuting differential operators pj''-^- 



p(i) 



N 

^ dW 
n {D, + 2tk), D,=p, + z—^, 

k=-M ^'ii 

0, [P}i\pj!;^] = 0, j,k = l,...,r. 



;8.9) 

^.10) 
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The A/'-fold supersymmetry is generated by 

in which ipj and are the annihilation and creation operators of the j-th fermion. 

We will not discuss the quasi-exact solvability of multiparticle rational A type Inozemtsev 
model, since its wavefunctions are not square integrable as seen in section |^. One could say 
that the A/'-fold supersymmetry is spontaneously broken in this case. 

Let us summarise that the quasi-exact solvability of quantum Inozemtsev models dis- 
cussed in sections H |^ and | is a consequence of the exact solvability of the quantum 
Calogero-Moser models and the quasi-exact solvability of the added single particle like in- 
teractions. 



9 Comments and discussion 

It should be stressed that the present method for showing quasi-exact solvability of single 
particle systems developed in section |^ to |^ does not depend on any existing methods or 
criteria for QES, see for example, As shown in section ^ it also gives the known exact 
solutions when the QES system tends to the harmonic oscillator. 

An interesting question along the line of arguments in this paper is the 'hierarchy problem' 
of quasi-exactly solvable systems, as in the completely integrable systems. For example, the 
Inozemtsev models have higher conserved quantities Tr(L'^) obtained from the Lax pairs in 
Appendix A. They define new classical and quantum Hamiltonian systems. Can the quantum 
version of the higher members of the hierarchy be deformed to be quasi-exactly solvable? 

in So far, elliptic Calogero-Moser models defied various attempts to construct quantum 
theory based on a Hilbert space, although existence of mutually commuting operators are 
known for the A type models |3l|]. In analogy with the present arguments, it is quite natural 



to expect the quantum elliptic Calogero-Moser models to be quasi-exactly solvable [^, ^ 
rather than exactly integrable. 
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Appendix A: Lax pairs for classical Inozemtsev models 

Here we present the Lax pairs for classical Inozemtsev models in the same notation as is used 
in the main text as evidence for their classical integrability. As mentioned in Introduction 
only certain subset of classical Inozemtsev models can be made quasi-exactly solvable. We 
focus on the supersymmetrisable Inozemtsev models for simplicity of presentation. For the 
full content of Lax pairs of classical inozemtsev models we refer to the original paper by 
Inozemtsev and Meshchryakov and for the universal Lax pairs for Calogero-Moser models 
in general, see 0, |TB[. 



The Lax pair consists of a pair of 2r x 2r (r rank) matrices L and M such that the 
canonical equations of motion can be expressed in a matrix form 

L = [L, M] 

and a sufficient number of classical conserved quantities can be obtained as as the trace of 
powers of L, Tr(L'^). The Hamiltonian is if oc Tr(L^). 

A.l BC type models 

The following Lax pair applies for the models presented in subsection p.5| for the rational as 
well as trigonometric models for proper choice of functions, x, u, etc. as listed below. The 
pair of matrices decomposes into diagonal and off-diagonal matrices: 

L = P + X, M = D + Y. (A.l) 

The diagonal matrices P and D are of the form 

r r 
P = "^PjiEjj - Ej+r,j+r), D = Dj{Ejj + Ej+r,j+r), (A. 2) 

j=i i=i 

in which Ej^^ is the usual matrix unit {Ej^k)im = ^ij^mk- The diagonal-free matrices X and 
Y have the form 

X = igu x{qj - qk)Ej^k + Wm Y1 x{qj + qk)Ej^k+r 
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+ 'i'9M ~ Qk)Ej+r,k + igM + (lk)Ej+r,k+r 

j¥=k 3+k 

+ 2i 5^ iy{qj)Ejj+r - 2i iy{qj)Ej+r,j, (A.3) 

i j 

y ^ y(^j ~ (ik)Ej,k + igu y(^j + (ik)Ej,k+r 

j¥=k j+k 

+ ^9m Y y(~1j ~ 1k)Ej+r,k + WM Y yi-(Ij + qk)Ej+r,k+r 
j¥=k 3+k 

+ ^ Y '''i(lj)Ej,j+r + i Y ^\<l3)E3+r,r (A.4) 

3 3 

The diagonal elements of D are given by 

r r 

Dj = -igM Y i^i^ij - ^k) + z{qj + qk)) - i Y ^(^i)- i^-^) 

k+3 J=l 

Some functions are related to each other: 

y{u) — dx{u)/du^ z{u) — x{uY + const., t{u) — 2x{2u)v{u) + const. (A. 6) 

The rational and trigonometric models correspond to the following choice of functions: 

1. Rational model, 

xiu) — — , ziu) — — , viu) — —iau^ + hu) + — , (A. 7) 

where a, 6 and 5^5 are real coupling constants. 

2. Trigonometric model, 

xiu) = cotw, ziu) = — 9 — , viu) = asin2ti — h occotti, (A. 8) 

where a, 6 and gs are real coupling constants. 
The functions x and correspond to those appearing in dW/ dq. 

A. 2 A type models 

The Lax pair is again a pair of 2r x 2r (r — 1 is the rank) matrices, with the decomposition 

L^P + X, M^D + Y, 
in which P and D are diagonal matrices of the form 

r r 

P = Y Pji^j,J ~ Ej+r,3+r), E) =Y Ej{Ejj + Ej^r,3+r) 
3=1 3=1 
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and X and Y are diagonal-free matrices of the form 



+ 2 E <^j)Ejj+r + 2 E <-Qj)Ej+r,j, (A.9) 
i j 

y = yiQj - qk)Ej,k + i9j2 yi^lj + C[k)Ej+r,k+r 

j^k i+k 
+ E + E «^'(-g,)^,+r,i. (A.IO) 

The diagonal elements of D are given by 

r 

= -ig E - Ik)- (A.ii) 

The rational and trigonometric models correspond to the following choice of functions: 

1. Rational model, 

x(u) = — , z(u) = — , k{u) = av? + 6m, (A. 12) 

u 

where a and b are arbitrary real coupling constants. 

2. Trigonometric model, 

x(u) = , z(u) = — 2 — , K,{u) = asin2M, (A. 13) 

sinw sin u 

where a is the arbitrary real coupling constant. 

Appendix B: Lower triangularity 

Here we show the details of the argument that Calogero-Moser part of the similarity trans- 
formed Hamiltonian (|7.10|) is lower triangular in the basis (|7.13|) . The triangularity of the 



trigonometric Calogero-Moser model is proved universally in ||T6[ by using the Coxeter (Weyl) 
invariant basis: 



in which A is a dominant weight 



A = E^^i-^i' ^ '■ fundamental weight (B.2) 
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and W{X) is the orbit of A by the action of the Weyl group: 



W{X) = {fie A(A)| fi = g{X), e G^}, A(A) : weight lattice of A. (B.3) 
The above (f)\{q) is Coxeter invariant. The set of functions {0a} has an order >: 

|A|2>|Af (pxxpy. (B.4) 

For the BCr root system, the set of weights W{\) is symmetric: 

/X G W{X) ^ G W{X). (B.5) 
Thus the Coxeter invariant basis ( p.l|) for BC type root system can be rewritten as 

<P'xi(l)= E cos(2/i-g). (B.6) 

All the fundamental weights, except for the spinor weights of are integral. That is, so long 
as the above dominant weight A contains the fundamental spinor weights in even multiples, 
all the /i ■ g in ( p.l|) have the form 

r 

fx-q = J2kjqj, kjEZ^. 



Therefore the basis ( |B.6| ) has the form ( [7.13| ) used in section |^. Moreover, after the applica- 



tion of H on the basis functions ( |7.13| ), those appearing as lower terms have also the same 
property of having even number of fundamental spinor weights and can be expressed in the 
same form ( |7.13D . This completes the proof of the lower triangularity of the Hamiltonian 
( |7.1(J| ) in the space ( |7.13| ). 
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